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Using scanning tunneling spectroscopy (STS) at 5 K in B-fields up to 7T, we investigate the local 
density of states of a two-dimensional electron system (2DES) created by Cs adsorption on p-type 
InSb(llO). The 2DES, which in contrast to previous STS studies exhibits a 2D Fermi level, shows 
standing waves at B = OT with corrugations decreasing with energy and with wave numbers in 
accordance with theory. In magnetic field percolating drift states are observed within the disorder 
broadened Landau levels. Due to the large electric field perpendicular to the surface, a beating 
pattern of the Landau levels is found and explained quantitatively by Rashba spin splitting within 
the lowest 2DES subband. The Rashba splitting does not contribute significantly to the standing 
wave patterns in accordance with theory. 

PACS numbers: 75.70.Tj, 73.21. Fg, 68.37.Ef, 71.20.Nr 



I. INTRODUCTION 

The spin-orbit coupling in semiconductors provides the 
opportunity to manipulate spins by electric fields^ which 
might become a central handle in spintronics^^i^ The cor- 
responding spin-splitting is called the Rashba effect j^*^ 
It has been probed in Ill-V-semiconductors, e.g., by the 
beating pattern of Shubnikov-de Haas oscillations-'^ or 
by the analysis of weak antilocalization.^^^ The Rashba 
effect has also been probed on metal surfaces by angular 
resolved photoelectron spectroscopyiSrJ^ as well as by 
scanning tunneling spectroscopy (STS)^i^ii^ In the lat- 
ter case, either the enhancement of the density of states 
close to the onset of the two-dimensional surface bandit 
or the suppression of quasiparticle interference caused by 
a missing spin-Umklapp scattering^^ has been used to de- 
duce the Rashba effect indirectly. Theoretical considera- 
tions reveal that the Rashba splitting is not directly visi- 
ble in quasiparticle interference patterns probed by STS, 
if only single scattering is considered. Subtle changes 
of the quasiparticle interference appear, if multiple scat- 
tering becomes relevant. ^- 

So far, Rashba spin splitting within semiconductors 
has not been probed in real space. Here, we use the 
surface doping effect, which leads to a two-dimensional 
electron system (2DES) directly at the surface of a 
low-gap III-V semiconductor to be probed by STS^i^Ti^i 
These 2DESs exhibit a quantum Hall effect down to 
B = 2T (filling factor 7) and a reasonable mobility of 
/i = 6000cm^/(Vs)j2^i2^ The low gap provides, in ad- 
dition, a relatively large Rashba splitting^^ as well as 
a strong Landau^^ and spin splitting.2^ This 2DES has 
been probed previously by STS revealing transitions from 
strong to weak localization,^ a spatially continuous wave 
pattern caused by multiple scattering, drift states in 
magnetic field^^ and the local density of states across 
quantum Hall transitions.^^ 

Within this paper, we describe STS measurements of 
a 2DES, which is induced by a minute amount of Cs 
(1.5 % of a monolayer) on the strongly p-doped InSb(llO) 



surface (A^a — 10^^ m~^). The strong doping results 
in a strong electric field (> lO^V/m) within the 2DES 
and correspondingly leads to a large Rashba coefficient 
{a IQ-i^eVm). We demonstrate that this 2DES still 
exhibits spatially continuous wave patterns at 5 = OT 
with preferential wave vectors in agreement with k • p- 
theory as well as strong corrugation (50% at £^f), very 
similar to the 2DES prepared on n-type InAs(l 10)^2^ 
Also Landau levels and drift states are observed by STS 
m a 5-field.2^ The disorder given by the acceptors pro- 
hibits the observation of spin splitting within the spa- 
tially averaged density of states, but spin splitting is 
observed in the local density of states {LDOS{x,y)). 
Importantly, the Rashba effect leads to a pronounced 
beating of the Landau level intensity measured by STS, 
which is quantitatively reproduced by calculations using 
a = 7 X 10~^^eVm. This experimental value is very 
close to the expected value of a = 9-11 x 10~^^eVm 
deduced from the known surface band bending. Thus, 
we demonstrate that the Rashba parameter can be de- 
termined down to the nm scale. 



II. EXPERIMENT 

Our home-built scanning tunneling microscope (STM) 
operates within an ultra-high vacuum (UHV) insert of a 
helium-4 bath cryostat at a base temperature of 5 Kelvin 
and with a magnetic field up to 7 Tesla perpendicular to 
the sample surface. For tip exchange and sample trans- 
fer, the STM can be lifted out of the cryostat into a 
standard UHV chamber without breaking the vacuum. 
A cesium dispenser (SAES Getters) mounted into this 
UHV chamber can dispense controlled amounts of ce- 
sium onto samples in the cold microscope. The STM tip 
was etched outside of the vacuum system from a tung- 
sten wire and prepared within the STM by field emis- 
sion and consecutive voltage pulses on a W(llO) crys- 
tal. As a sample we used a Ge-doped InSb crystal with 
a hole concentration of 1-2 x 10^^ m~^ as determined 
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by Hall measurements. In order to prepare the 2DES, 
the crystal was glued to a molybdenum sample holder 
with electrically conductive adhesive pointing with the 
(110) surface upwards. It was cleaved at room temper- 
ature at a background pressure of p = 2 x 10~^^mbar 
resulting in a clean (110) surface. The sample was then 
transferred into the pre-cooled STM (T < 60 K) and ce- 
sium was evaporated onto the crystal surface. Without 
cooling, the adsorbed cesium atoms would quickly form 
chain structures at room temperature^^ breaking the ho- 
mogeneous distribution. After transferring the micro- 
scope into the cryostat the sample could be probed for 
weeks without any noticeable change in adsorbate den- 
sity. The cesium atoms appear as white dots in the STM 
image of Fig.[2fa), which shows a homogeneous Cs distri- 
bution amounting to about 1.5% adsorbates per surface 
unit cell. STS data {dl/dV) were acquired by lock-in- 
technique applying a modulation voltage with amplitude 
a/2 X Vmod to the sample. Prior to measurement, the tip- 
surface distance is stabilized at voltage Fgtab and current 
^stab- Then the feedback is turned off and the voltage is 
ramped to the measurement voltage Vg. 



III. RESULTS AND DISCUSSION 

It is well known that adsorbates like alkali metals can 
induce band bending at semiconductor surfaces already 
at low adsorbate densities. A strong downwards band 
bending on p-type semiconductors induces a so-called 
inversion- layer, which confines a 2DES directly below 
the surface. In contrast to 2DESs buried deeply in 
heterostructures,^^ this 2DES is accessible by STS, a 
technique, which measures the LDOS{x^y) by the dif- 
ferential conductivity dl / dV {x ^ y) .^^ 

The strength of the adsorbate induced band bend- 
ing mainly depends on the materials combination and 
the adsorbate coverage. ^^^^ Cs on n-InSb(llO) induces 
a maximum band bending of Fbb = 290 meV as deter- 
mined by photoelectron spectroscopy,^^ where Vbb de- 
scribes the energy difference between the Fermi level 
and the onset of the conduction band directly at the sur- 
face EcBM' Therefore, the maximum band bending for 
Cs on p-InSb(llO) is Vbb + £^gap with the energy gap 
Egap = 0.235 eV^^ of InSb. The actual shape of the con- 
fining potential and the resulting subband energies of the 
2DES further depend on the dopant concentration of the 
semiconductor. The steepness of the band bending in- 
creases with the acceptor concentration Na- A large A^a 
results in a steep confining potential and, thus, leads to 
large subband energies Ei, i = 0,1,2,... measured with 
respect to £^cbm- Consequently, one gets a low electron 
concentration within the inversion-layer N^. 

Assuming the subband energy levels Ei to be known, 
the electron concentration in the inversion- layer is 
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FIG. 1. Spatially averaged dl/dV curves measured at B = 
T and B = 7 T as indicated, V^tab = 300 mV, /stab = 200 pA, 
Vmod = 1 mVrms. 12 X 12 spcctia were taken over an area of 
300 nm x 300 nm. The subbands Eq and Ei are marked at 
— 60mV and 110 mV, respectively. 



given by 
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where D{E) denotes the two-dimensional density of 
states including spin degeneracy and m*{E) is the ef- 
fective electron mass within the conduction band, which, 
in general, is energy dependent due to non-parabolicity. 
An estimate for the so-called momentum effective mass 
or density of states effective mass for a 2DES defined 

oo36,37 
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can be given within a triangular potential well 
approximation^^ using a k-p perturbation approach:— 
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Here, E splits up into the in-plane energy E\\ = E — Ei 
and the energy offset of the subband Ei with respect to 
EcBM^ ^0 denotes the bulk effective electron mass at 
^CBM and £^gap is the gap energy. Strictly speaking, 
this approach is limited to energies small compared to 
the gap energy. But although this is not the case for the 
2DES induced by Cs on InSb, Eq. dH can stih be used 
to get a reasonable estimate of the effective mass. 

Figure [1] (lower curve) shows the dI/dV{V) spectrum 
of Cs covered InSb(llO) at 5 = OT, which is averaged 
from 144 curves recorded on a regular grid covering an 
area of (300 nm)^. The first and second subband can 
be identified as steps in the dl/dV signal, which due to 
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FIG. 2. (a) STM constant-current image of InSb(llO) covered with 1.5 % of a monolayer Cs; Cs adsorbates appear as white 
dots; the size of the inset is indicated by the white rectangle; 200 nm x 200 nm, V = 300 mV, / = 100 pA; (b)-(j) dI/dV{x,y) 
images of the same surface area as shown in (a) recorded at F = —60 mV to V = 20 mV as marked in the insets; V^stab = 300 mV, 
^stab = 100 pA, Vmod — 1 HiVrms, pixcl rcsolutiou: 2 um; insets show Fourier transformations of the real-space images displaying 
the electron wave vectors contributing to the LDOS{x, y). 



averaging represents the density of states (DOS). The 
applied sample voltage corresponds to the electron 
energy E relative to the Fermi energy Ey- Obviously, 
only the first subband is occupied by electrons, since the 
second subband is located above E^ iV^ = OmV). As- 
suming reasonably that the band bending induced by Cs 
is complete at 1.5% coverage, i.e. Vbb = 290 meV,— 
the subband energies Ei relative to £^cbm are deduced 
to be £^0 = 230 meV and Ei = 400 meV. Using Eq. (g]) 
with tuq = 0.0135 x mo^^ (mo: free electron mass), the 
effective electron mass in the first subband increases from 
0.022 X mo at the subband edge {Eq) towards 0.029 x mo 
at the Fermi level (^f)- The effective mass at the onset 
of the second subband {Ei) would be 0.042 x mo- Using 
further Eqs. ^ and (|2]), the electron concentration of the 
2DES can be deduced to be = 6.5 x 10^^ m-^ No- 
tice, that the dl/dV curve of Fig. [T]( lower curve) does not 
show any states of a tip induced quantum dot,^^ which 
has been achieved by careful preparation of a tip with 
adaptive work function to the sample surface. 

The dI/dV{x^y) (LDOS) images probed at energies 
within the first subband are shown in Fig. [2fb)-(j). One 
observes the evolution of standing wave patterns exhibit- 
ing decreasing wave length with increasing energy. These 
patterns are very similar to the ones observed within a 
2DES induced by Fe on n-type InAs(llO).^^ Notice that 
the small black dots appearing in all dl/dV images at 
the same position are the Cs atoms, which appear dark 
due to the larger distance of the tip with respect to the 
InSb surface. The insets show the Fast Fourier Transfor- 
mations (FFTs) of the real-space images displaying the 
electron wave vectors k contributing to the real-space 
pattern of the LDOS with k = 0/nm being in the cen- 
ter. In the ideal case of a 2DES with negligible poten- 




FIG. 3. (a), (b) dI/dV{x,y) images recorded at the Fermi 
energy {Vs = OmV) and at B = OT (a) and B = 7T 
(b); 300 nm x 300 nm, Vst^h = 300 mV, /stab = 200 pA, 
Vmod = ImVrms, plxcl resolutlou: 1 nm; insets show Fourier 
transformations of real space data. 

tial disorder, the EFT would show a ring, growing in di- 
ameter with the non-parabolic energy dispersion relation 
E{k). The contribution of larger k values with increas- 
ing energy can indeed be deduced from the growing disc 
in the FFTs at low energy, which develops into a ring 
structure at about Vs = mV still increasing in diameter 
with energy. The contrast of the ring can be improved 
by recording larger images with a higher spatial resolu- 
tion as shown in Fig. [3l^a) for the wave pattern at Ep. 
A clear ring structure appears in the EFT image, which 
is difficult to observe in Fig. Efh). The transition from 
a disk-like appearance towards a ring-structure with re- 
maining intensity in the center of the ring has been ob- 
served previously for a disordered 2DES.^^ There, it has 
been reproduced by Hartree calculations taking into ac- 
count the potential disorder of the 2DES, produced by 
the charged dopants. For strongly disordered systems, 
one observes only the disk increasing in diameter with 
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FIG. 4. The electron energy E relative to the Fermi energy Ef 
is plotted against the k value parallel to the 2DES layer. The 
dots represent the dominant k values extracted from Fig. [2] 
(e)-(j). The upper line shows a parabolic dispersion with the 
estimated effective mass of m* = 0.022 x mo at the subband 
edge. The lower curve is a numerical solution of Eq. Q and 
Eq. (|4|). The deviation from the parabolic dispersion marks 
the non-parabolicity of the 2D conduction band dispersion 
of InSb, the deviation from the numerical approximation is 
caused by the influence of higher order terms in k • p theory. 



energy.^^ Thus, the FFTs, not being a perfect ring struc- 
tures, indicate the wave function mixing by the spatially 
fluctuating electric potential due to charged dopants. In 
accordance with expectations, the wave function mixing 
gets reduced with increasing energy. 

Figure [H shows the dominant k values taken from the 
maxima in radial line scans of the FFT images in Fig. [2] 
and [3l They are displayed as dots in comparison with 
a theoretical upper and lower limit of the expected InSb 
dispersion. The upper limit is a parabolic dispersion re- 
sulting from the effective mass m* = 0.022 x mo at the 
onset of the 2DES Eq. It neglects the non-parabolicity 
within the 2DES. The lower curve is obtained by solving 
Eq. (j3|) and Eq. (|3j) numerically, which overestimates the 
non-parabolicity at the high energies of the 2D electrons 
by neglecting higher order terms of k-p theory.^^^^^ A full 
numerical treatment of the k • p Kane model^^ is beyond 
the scope of this paper. Nevertheless, the experimental 
data points are found in between the two limits evidenc- 
ing that the wave patterns of Fig. [2] are indeed caused by 
the electrons of the 2DES. 

Interestingly, the wave pattern intensity is rather con- 
tinuously distributed over the area of Fig. [2] and [3l^a). 
This can be most clearly seen in Fig. [3fa) and is in con- 
tradiction to expectations from single point scattering, 
where a reduction of intensity / with distance r from the 
scatterer according to /(r) ex is expected. One rea- 
son might be that this discrepancy is influenced by multi- 
ple scattering paths to the standing wave pattern, which 
are known to be considerably more important in 2D than 
m 3D 4^ Another reason might be the large number of 
scatterers within the image area of Fig. [2j This number 
can be estimated with the help of the depth of the 2DES 
as displayed in Fig. [8] (9nm), the acceptor concentra- 
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FIG. 5. (a) Corrugations K of the electron wave patterns 

within the first subband; (1) K := (Cmean - Cmin)/Cmean, 
(2) K (2 X Cstd)/Cmean with Cstd, Cmean and Cmin bclug 

standard deviation, mean and minimum in the intensity dis- 
tribution of the dl/dV images, (b) Histogram of Fig. EJJh) 
indicating the derived values. Cmin marks the dl/dV offset 
from a linear fit of the left slope (dotted line) of the histogram 
which amounts to cutting off the lowest 2% of dl/dV values. 

C'mean and Cstd 

are calculated numerically. 



tion (10^^ m ^) and the image size ((200 nm)^). This re- 
sults in 360 acceptors scattering the electron waves within 
Fig. El 

The corrugation K of the wave patterns, which is 
caused by the strength of the scattering and the phase 
coherence length, is displayed in Fig. [5fa). It is calcu- 
lated by two different methods using the histograms of 
dl/dV values within an image as shown in Fig. [5fb). 
The first method determines the difference between the 
mean value Cmean and the minimum value Cmin and di- 
vides it by Cmean- K = (Cmean ~ C^min) / Cmean' The 

second method instead uses twice the standard devi- 
ation of the histogram Cstd and divides it by Cmean- 

= 2 X Cstd/C'mean- Both approaches give similar results 
except at low energy as visible in Fig. [5fa). In previous 
STS studies, method (1) has been used revealing that a 
2DES {K = 60 %) is prone to a much larger i^- value than 
a three-dimensional system {K = 3%), even if the po- 
tential disorder is quite similar. Moreover, a drop of K 
from 90 % to 50 % has been found in a strongly disordered 
2D system at the percolation transition of strongly local- 
ized states. Our i^-values are compatible with the pre- 
vious values obtained on 2DESs, but, other than in the 
previous publications, K becomes continuously smaller 
with increasing energy. This can be explained straight- 
forwardly by the increase of the scattering length with 
energy. However, the difference with respect to the pre- 
vious data is not completely clear. 

Next, we discuss the STS data obtained in a perpen- 
dicular magnetic field of B = 7T. The dl/dV image 
recorded at the same sample position as Fig. [3](a) us- 
ing exactly the same tunneling parameters is shown in 
Fig. [3](b). The pattern looks much more disordered and 
a dominant wave length is not visible anymore as ev- 
idenced by the FFT in the inset. In a perpendicular 
magnetic field. Landau quantization and spin splitting is 
expected for 2D electrons. In the effective mass approxi- 
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FIG. 6. (a)-(f) dl/dV images with 1.5 nm resolution at B = 
7T, 150nmxl50nm, V^tab = 300 mV, /stab = 200 pA, K,od = 
ImVrms; the corresponding sample voltages are marked in 
the local dl/dV curves shown in (g); (g) local dl/dV curves 
spatially averaged over the two small areas marked in the 
images (a)-(f). 



FIG. 7. (a)-(f) dl/dV images recorded at B = 7T, V^tab = 
300 mV, /stab — 200 pA, ymod — ImVrms, with different sam- 
ple voltages as marked in (g), 150 nm x 150 nm; the intensity 
marked by black lines highlights a drift state moving uphill 
with increasing energy; (g) spatially averaged dl/dV spec- 
trum originating from the 10^ curves covering the whole im- 
age area of (a)-(f). 



mation neglecting Rashba spin splitting, the Landau and 
spin levels are given by 

=E,^ hujcin + ± ^^*/iB5 (5) 

with subband index z, Landau level index n = 0,l,2,..., 
spin index ±, Bohr magneton fiB = eh/2mo, elementary 
charge e, cyclotron frequency 



eB_ 



(6) 



and effective Lande g-factor The spin splitting given 
by the third term in Eq. ^ is also energy dependent 
within a non-parabolic conduction band, i.e. g^{E). An 
approximation for g^{E) is given by the relation:^^- 



9*iE) 
9o 



m*{E)' 



(7) 



For InSb, one can use = —51 being the effective g- 
factor at the conduction band edge.^^ 

Within an electric potential landscape, the Landau 
states experience an additional electric field perpendic- 
ular to the 5-field, which leads to drift states, i.e. states 
covering closed equipotential lines of the disorder. ^^^^^ At 
the lowest possible energy, the first Landau level starts 



with states at potential pits representing localized elec- 
trons. With increasing energy, the equipotential lines 
representing the states encompass an increasing area un- 
til they percolate at the critical energy close to the center 
of the Landau level. There, the extended critical state 
of the Quantum Hall transition appears. At even higher 
energy, the states localize again by covering equipoten- 
tial lines circulating around potential maxima. This be- 
havior, taking place in each Landau and spin level, has 
recently been observed directly by STS at T = 0.3K.^°- 
A similar percolation of states in 5-field has also been 
observed on graphite surfaces. 

For our system, the transition is shown in Figs. [6] and 
[3 concentrating on the lowest Landau level. The im- 
ages (a)-(f) in Fig. [6] display dl/dV images at 5 = 7T 
at increasing sample voltage. The circle drawn in each 
image marks a potential pit, since it shows dl/dV in- 
tensity already at the lowest energies. The square, in 
contrast, marks an area, where the potential is relatively 
high, i.e. the dl/dV intensity appears only at rather high 
energy. The curves taken at these positions are shown in 
Fig. [6fg) and indeed exhibit a pair of peaks shifted by 
approximately 20meV with respect to each other. One 
can identify the doublet as the first spin split Landau 
level, which would result in |^*| = 32 for the spin split- 
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ting and m* = 0.020 x mo for the distance to the next 
Landau level observed within the pit. This is in excellent 
agreement with the values \g* \ = 31 and m* = 0.022 x mo 
calculated for the subband edge Eq from first-order k • p- 
theory (Eqs.|4]and ([Tj)). By comparing the two curves of 
Fig. El^g), we deduce that our 2DES exhibits a disorder 
potential with a potential fluctuation of about 20meV 
and, by regarding the distance of square and circle in 
Fig. [6fa)-(e), we deduce a spatial length scale of the po- 
tential corrugation of about 30 nm. 

Since the potential fluctuation is larger than the spin 
splitting, it is difficult to observe the appearance and dis- 
appearance of extended states. However, a first extended 
state is visible in Fig. El^c), where the states developing 
from the potential pits touch and form connecting paths 
from one side of the image to the other. An upward move- 
ment of a drift state onto a potential hill is highlighted 
in Fig. [3 The black lines roughly mark the dl/dV in- 
tensity of a closed structure, which at low energy sur- 
rounds an area of about (80 nm)^. With increasing en- 
ergy, the dark inner area shrinks continuously condensing 
towards a small spot in (f). The overall increase of the 
spatially averaged dl/dV signal in (g) can be attributed 
to the drift states of the next Landau level, which start 
to develop in the lower potential areas. The visibility 
is disturbed by the overlap with states from other lev- 
els. Therefore, Fig. [Tl^a) is taken from the lower spin 
level and Figs. [7fb)-(f) are taken from the higher spin 
level leading to reduced overlap with other states in both 
cases. Indeed, the energetic overlap of states from differ- 
ent Landau and spin levels at different potential energy 
completely removes the spin splitting from the spatially 
averaged dl/dV curve, i.e. the DOS, and weakens the 
visibility of the Landau level splitting significantly, as 
shown in Fig. [3g). 

Finally, we discuss the observation of the Rashba effect 
within our 2DES. As described in Eq. ([5]), constant values 
of m* and would lead to equally spaced Landau and 
spin levels at constant 5-field. Each level is broadened 
by the disorder potential, which is relatively large in our 
system because of the high acceptor concentration pro- 
hibiting the observation of spin splitting in the spatially 
averaged density of states. For non-parabolic conduc- 
tion bands, the Landau level splitting huJc{E) and the 
spin splitting |^*(^)|/iB^, both, decrease with energy. 
This fact, however, cannot explain the spatially averaged 
spectrum shown in Fig. [T] (upper curve), which exhibits a 
beating pattern of the Landau level intensity within the 
first subband. The proposed explanation for the beating 
is Rashba spin splitting,^ which gives a spin splitting for 
confined 2D electrons moving perpendicular to an elec- 
tric field even at 5 = T. The electric field can be exter- 
nally applied or is given by an asymmetric confinement 
potential. The Rashba spin splitting increases with wave 
number /c, leading to two dispersion curves with differ- 
ent effective mass. Including this effect into the effective 
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FIG. 8. Calculated band bending for the 2DES: conduction 
band minimum (CBM), valence band maximum (VBM) and 
subband energies Eq , Ei are plotted as a function of distance 
z from the surface. The assumed acceptor concentration and 
electron concentration of the 2DES are Na = 1 x 10^^ m""^ and 
Ns — 6.5 X 10^^ m~^, respectively, in accordance with experi- 
mental data. The electron distribution curves |^i(2;)|^ result- 
ing from a triangular well approximation are drawn for each 
subband. "^"^ The depletion length Zdepi = 29 nm and the elec- 
tric field visible within the first subband \E\ = 3.1 x 10^ V/m 
are additionally marked. 



mass Hamiltonian results in Landau levels given by^-- 

= Ei + TiLOe (n + a {5^ + 7'n) , (8) 
^ = a(2m*/K'uj,f'\ (9) 

(10) 



2 V 2mo 



with n = 0, 1, 2, . . . being the level index, a = 1 for n = 0, 
and cr = ±1 for n = 1, 2, 3, . . . being the spin index. The 
parameter a is called the Rashba parameter. Here it en- 
codes the interaction of the conduction band with the 
spin-orbit split valence band and can be deduced from 
k • p-theory.^^^^^ One can easily show that Eq. (|8]) gives 
the same energy levels as Eq. ([5]) for a = 0. It is obvi- 
ous from Eq. ([8|), that the g- factor induced spin splitting 
given by 5 is the same for all n, while the a dependent 
splitting encoded in 7 increases with the level index n. 
This eventually leads to a mixing of spin states from dif- 
ferent n in the density of states inducing a beating. 

The Rashba parameter a for the 2D conduction band 
electrons of InSb can be estimated in an eight band Kane 
model^^ with a confining potential V{z),^^ thereby de- 
scribing non-parabolicity as the interaction of the spin 
split conduction band with the nearest three spin split 
valence bands, of which one is separated by the spin- 
orbit splitting energy A = 0.8 eV. This results in an 
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FIG. 9. Spatially averaged dl /dV (DOS) spectra at different B-fields as marked (experiment) in comparison with simulations 
using Landau and spin energies from Eq. (|8]) with m* = 0.035 x mo, a — 7 x 10~^^ eVm, ^* = —21 and a Gaussian broadening 
with full width at half maximum (FWHM) as indicated for each B-field; the beating minima are marked by arrows. 



energy and z dependent Rashba parameter jsi 

P2 d / 1 

a[z,E) - — — \E-V{z)^Eg 



E-V{z)^Eg^Aj ' ' ^ 

3Eg^2A ' ^ ^ 

To calculate the spin splitting via Eq. (|8]), a{z^E) 
has to be convolved with the wave function of the low- 
est subband in z-direction ^o(^)- In order to obtain 
^0(^)5 the confining potential is calculated as shown in 
Fig. [H For this calculation, the maximum cesium in- 
duced band bending Vbb = 290 mV has been used to 



solve the Poisson equation. The screening charge is given 
by the bulk acceptor concentration A/'a, which is nega- 
tively charged by electrons from the adsorbates, and the 
estimated inversion- layer electron concentration = 
6.5 X 10^^ m~^.^^ We find, that A^s has a negligible influ- 
ence on the resulting confinement due to the dominating 
2D acceptor concentration Na x Zdepi = 3.0 x 10^^ m~^. 
Using the triangular well approximation with an infinite 
barrier at 2; = 0nm4^ and the effective mass given by 
Eq. dH), two subbands result within the confinement po- 
tential exhibiting energies in good agreement with the 
experimental values shown in Fig. [H More precisely, 
the calculated energies £^0 and Ei are by about 10 % 
larger than the experimental ones, which is not surpris- 
ing, because the simple model neglects penetration of 
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the electrons into the vacuum and into the valence band, 
which both would lower their energy. A more accurate 
calculation would require solving the Poisson equation 
and a multi-band Kane Hamiltonian self-consistently,^^ 
but this is beyond the scope of this paper. Calculating 
the band bending for the whole width of possible accep- 
tor concentrations Na = 1-2 x 10^^ m~^ reveals electric 
fields within the 2DES of |^| = 3.1-4.1 x 10^ V/m. E as 
well as the square of the wavefunctions result- 
ing from the triangular well approximation are drawn in 
Fig. [8] highlighting the nearly constant value of E across 

Finally, the effective Rashba parameter a for the first 
subband is determined convolving the result of Eq. ([TT]) 
with ^o(^) numerically: 

a={^o{z)\a{z,Eo)\Mz))z. (13) 

This results in a Rashba parameter a = 9-11 x 
10~^^ eV m being higher than the values observed in In As 
inversion- layers or heterostructures by transport mea- 
surements (3-4 X 10~^^ eVm).^^ It should be noted that 
the calculated Rashba parameter is an upper estimate, 
since the ignored effect of barrier penetration of the elec- 
tronic waves decreases the effective Rashba parameter. 
Furthermore, a is only the lowest order of an inversion 
asymmetry induced spin splitting and it is known that 
higher orders lead to a reduced effect. Indeed, the ob- 
served beating patterns displayed for different 5-fields in 
Fig. [9] can be reproduced nicely by using Eq. (j8j) with a 
slightly reduced Rashba parameter ofa = 7x 10~^^ eVm. 
For the sake of simplicity, we assumed a constant effec- 
tive mass and ^*-factor taken as the average value within 
the first subband (m* = 0.035 x mo, g* = -21) and we 
fit a Gaussian broadening to each spin and Landau level 
accounting for the potential disorder in order to obtain 
the calculated curves in Fig. [9l Note that the fit param- 
eter FWHM = 13-18 meV given in Fig. [9] nicely agrees 
with the potential disorder deduced independently from 
Fig. El 

Although we could not reproduce all details of the mea- 
sured spectra, it is obvious that the minima of the beat- 
ing (nodes), marked by arrows in Fig. [9l are in excellent 
agreement with the experiment. The calculation also re- 
produces the increase of beat frequency with decreasing 
5-field. The decreased FWHM found for lower B re- 
flects the fact that the lateral extension d of the drift 
states scales according to d ex Thus, drift 

states become more insensitive to the steepest parts of 
the fluctuating disorder potential at lower B. We checked 
carefully, that the chosen FWHM does not influence the 
position of the beating nodes. 

The Rashba spin splitting is, of course, also present 
at 5 = OT and would split an approximated parabolic 
energy dispersion E{k) into two branches separated by 



if higher order spin splitting terms are neglected. 
This leads to a spin splitting of about 30meV at E^ 
and two /cp values of about k- = 2.7 x 10^ m~^ and 
/c+ = 2.1 X 10^ m~^. However, these two /cp -values are not 
visible in the EFT of Fig. [3] in accordance with theory. 
Slight changes within the complex LDOS pattern might 
appear, if multiple scattering is involved, but these 
changes could only be pinpointed by detailed compari- 
son with LDOS calculations including the details of the 
disorder potential. A more detailed theoretical consid- 
eration of the consequences of disorder on the Rashba- 
split LDOS based on previous calculations^'''^^ might be 
an interesting base for such experiments. Notice that 
the width of the ring in the EFT of Fig. [3fa) is as large 
as the difference of /c+ and k- at Ep showing that the 
Rashba effect and the disorder within our sample are of 
the same strength. We also checked, if a direct observa- 
tion of zero field splitting is possible within the density 
of states as has been recently proposed by Ref. [5l|, but 
without success. 



IV. CONCLUSION 

In summary, we have successfully prepared a 2DES 
close to the surface of a highly doped p-type InSb crys- 
tal exhibiting one occupied subband only. Using scan- 
ning tunnehng spectroscopy, the evolution of standing 
waves in the 2DES has been imaged being dominated 
by wave numbers explained by a non-parabolic disper- 
sion relation, but exhibiting strong wave function mix- 
ing and large corrugation due to disorder. In magnetic 
field. Landau and spin levels are observed locally, but 
washed out in the density of states due to the disor- 
der with amplitudes of about 20 meV. Percolating drift 
states are observed within the potential landscape. The 
density of states shows an irregular Landau level pat- 
tern dominated by beating which could be attributed to 
Rashba spin splitting caused by the asymmetry of the 
confining potential. The deduced Rashba parameter of 
a = 7 X 10~^^ eVm is relatively large and very close to 
the value of a = 9-11 x 10~^^ eVm estimated by k • p- 
theory. This shows that Rashba parameters can also be 
determined by a local probe, which provides spatial reso- 
lution down to the atomic scale. This opportunity might 
trigger novel types of experiments within the near future. 
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